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Abstract 

The essence of romance is mystery. In this talk, given in honor of the 60th 
birthday of Michio Jimbo, I will explore the meaning of this for the Ising model 
beginning in 1946 with Bruria Kaufman and Willis Lamb, continuing with the 
wedding by Jimbo and Miwa in 1980 of the Ising model with the Painleve VI 
equation which had been first discovered by Picard in 1889. I will conclude with 
the current fascination of the magnetic susceptibility and explore some of the 
mysteries still outstanding. 

1 Introduction 

A search of Google books reveals that the observation 

The essence of romance is mystery 

has been made by many authors in many different ways and in many different 
contexts ranging from the literary to the scientific. But in all contexts romance 
betokens fascination and the Ising model has fascinated many people, including 
myself, for many decades and in spite of many breakthroughs and moments of 
understanding the mystery continues to this day. In this talk I will present some 
of the milestones of this romance. 

2 Kaufman and Lamb 

In his talk "The Ising model in two dimensions" [I] presented at the fifth Battelle 
Colloquium on Materials Science, held in Geneva and Gstaad, Switzerland, 
September 7-12, 1970, Lars Onsager wrote, following a discussion of his famous 
1944 computation of the free energy [2] and a sketch of his 1945 proof of his 
conjectured spectrum of the transfer matrix, 



"Before long, however, Bruria Kaufman had developed a much better strat- 
egy 

At Columbia University she first asked Willis E. Lamb to direct her work 
on order-disorder problems; but he was much too heavily engaged in an exper- 
imental effort, and I was asked to assume the responsibility. Unable to talk 
her out of the idea I suggested that she explore • • • By the summer of 1946 she 
had a beautifully compact computation of the partition function, bypassing all 
tedious detail. 
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By itself that was only a more elegant derivation of an old result but the 
approach looked powerful enough to produce a few more new ones. Very well, 
how about correlations?" 

The history of the Ising model from that time forth has been the study of 
these correlations. 

But the deeper meaning of this passage from Onsager's paper completely 
escaped me until many years later Rodney Baxter wrote to me concerning a 
typescript [3] that had been given to him which is certainly a draft of Onsager 
and Kaufman's calculation of the spontaneous magnetization of the Ising model. 
Why in the world would Kaufman, who was creating pioneering mathematics, 
ask Lamb, an experimental physicist, to supervise her research? This question 
was brought into sharp focus when Baxter told me that he was going to contact 
her about the authorship of the typescript. She was then living in Tucson, 
Arizona with her husband, Willis Lamb. 

So this is the first romance concerned with the Ising model. Both Bruria 
and Willis were married to other people in 1946 when Bruria asked Willis to 
be her research supervisor and he turned her down. But decades later, when 
Kaufman's husband died in 1992, Lamb invited her to Tucson as a Visiting 
Scholar at the University of Arizona where he was a professor. In 1996, after 
his wife died, Willis and Bruria were married. 

3 Correlations and form factors 

The great understanding of Kaufman was that the Ising partition function could 
be written by use of fermionic methods as the sum of four Pfaffians [4] and that 
this fermionic method is powerful enough to write all correlation functions of 
the Ising model as determinants [S]. 

The Ising model is a system of "spins" <jj r k at row j and column k of a 
square lattice which take on the values tjj.k — ±1 and interact with their nearest 
neighbors with the interaction energy 

V L h 

£ = - E { Eh < J 3,kV J +i,k + E v a^ k <r j+hk }. (1) 

j=-L v k=-L h 

The correlation functions studied by Kaufman and Onsager are defined as 

(co.oCM.iv) = lim Zj} Lh V" ao,o&M,Ne~ £/kBT (2) 

L v ,L h — »oo ' — * 
ff=±l 

where T is the temperature, kg is Boltzmann's constant, 

= E e ~ £/kBT ( 3 ) 

CT=±1 

is the partition function and the sum X)o-=±i ' s over ai l values of the variables 

(Tj,k- 
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The discovery of Kaufman and Onsager [5] is that the row and diagonal 
correlations can be written as a sum of two determinants. These are further 
simplified by Montroll, Potts and Ward [5] to a single determinant. The diagonal 
(t T o,oCAf,A') and the row correlations (oo.oCo.Ar) can both be written as N x N 
Toeplitz determinants 



D 



N 



«0 
«1 



a_i 
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where 



with 



ttJV-l OJV-2 



a-N+i 

a -N+2 
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(l- ai e w )(l-a 2 e- w y 1/2 



(1 - aie- w ){l - a 2 e l6 )_ 
For (a ,ocr NtN ) 

ai = 0, a 2 = (sinh2E v /k B T sinh2E h / ksT)- 1 
and for (cro,oO'o,jv) 

ai = e - 2B " / ' £BT tanh J B'V/c B T, a 2 = e" 2£; " /fcBT coth^/fesT 



(4) 

(5) 
(6) 

(7) 
(8) 



and the square roots are defined to be positive at 6 = tt. These determinants 
are very efficient for the calculation of the correlations when TV is small. 

However, when iV is large the determinental representation ((4] is not an 
efficient method of calculation and a different representation must be found. 

The first step in finding this new representation is the computation of the 
limiting value as N — > oo 



lim (cro,o0-o,jv) 

N— too 



lim (<7 O ,o0'jv,iv) = ( 1 _ *) 1/4 

N—>oo 



with 



t = (smh2E v /k B Tsiiih2E h /kBTy 



(9) 



(10) 



which is valid for < t < 1. For t > 1 the limit vanishes. The value of T 
for which t = 1 is called the critical temperature T c . It is the evaluation of 
this limit for (ero.oCT/v.Ar) which is accomplished by Kaufman and Onsager in the 
manuscript recently published by Baxter 

The next step in the evaluation of the long distance behavior of the corre- 
lations was made in 1966 by Wu ;7] who computed the first correction /q fa to 
((9]) as N — > oo for (oo.oco.iv) f° r T < T c as a two-dimensional integral and the 
leading behavior /q 1 ^ as TV — > oo of (cto,oCo,at) for T > T c as a one-dimensional 
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integral. These are the first terms in what is now called the form factor expan- 
sion of the correlation functions, which for general M, N is written for T < T c 

as 

oo 

(<to, <tm,n) = (1 - t) 1/4 {l + £ (11) 

71=1 

and for T > T c as 

oo 

( ( To,oa M ^) = (l-t) 1/4 E/Mr^ 1) ' ( 12 ) 

where for T > T c we use the definition 

i = (sinh2£7fc B Tsinh2£;'7fc i 3T) 2 . (13) 

The derivation of the complete expansions ([111) and (|12l) has its own inter- 
esting story. In 1976 Wu, McCoy, Tracy and Barouch [5] derived an expansion 
valid for all N of the correlations in the form for T < T c of 

oo 

2n) 



(*o,octm,n) = (1 - <) 1/4 exp ^ jg£ (14) 



n=0 

and for T > T r 



(<to,o*m,n) = (1 - t) 1/4 £ exp E Ffift (15) 

n=0 n=0 

where F^ n ^ and F^ n ^ are 4n dimensional integrals and GK 1 ' are 4n + 2 
dimensional integrals. For all three functions half of the integrals may be exe- 
cuted by closing a contour integral on a pole. The forms (fT4"]) and (|15p of the 
correlation functions are called the exponential forms. 

The form factor expansions (fTTj) and (|12p are obtained from the exponential 
forms (fT^|) and (fl"5|) by expanding the exponentials. For a few low values of n 
this was done in [5] in connection with the study of the magnetic susceptibility 
but the general results for the f^ N were not given by Nickel [5] and [TU] until 
1999 and 2000. 

A curious feature of the derivation given in [5] of (fl"4"j) and (fT5")) is that the 
method of [7] developed for the row correlation (<to,o0o,n) is not used; instead 
the method used by Cheng and Wu [11] in the study of the leading terms of large 
separation behavior of the general correlation (ao t o(TM,N} is used. The original 
method [7] of Wu as applied to the correlations (cro 1 o0'o,iv) and (ao,oCN,N) was 
extended to all orders in 2007 by Lyberg and McCoy [12] . The results in [T2] 
for the diagonal form factors f^ N (t) are for T < T c 

t^ N+n) f 1 n rr ( (i-^-W-i) V /2 



4 



n n ( -i _ frr l, t^~ ) n (^-l-^fc-i) 2 ^-^) 2 , 

(16) 



y l - tx 2 k-lX2j , , 
l<3<n l<k<n x J ' !<3<k<n 



-l _m-i/2 



and for T > T c 

,(2n+l),,s _ 

t {n+l/2)N+n{n+l) i-l 2n + 1 n + 1 

n u n+1Vjr 2n+i / n n x 27-i [(! - - !)] 

' V 1/0 k=l j=l 

n , 1 

n^[ ( i- toJj ) (I --i)r ii n t^-^i 

j=l l<j<n+l l<fe<n v K 

Yl (x 2 j-l - X 2 k-l) 2 Yl (X 2 j-X 2 k) 2 - (17) 

l<j<fc<ri+l \<j<k<n 

A closely related form for the row form factor /i 7 ^ is also obtained in [TJ] . The 
results (fT6|) and (fT7)) have the startling feature that in the diagonal case the 
ZiViV ^° no * manifestly reduce term by term to the corresponding functions 
obtained from |8j. The reconciliation of these two forms is one of the present 
mysteries of the Ising model. 

These diagonal form factor integrals, which on the surface may appear to be 
indigestible, have proven to have many very special properties. 

1). All the integrals in (fTBT) and (|T7| reduce at t = to a product of two 
special cases of the celebrated Selberg integral [T3] 



/l pi n 

■••/ n*? -1 ^-**)"" 1 n \t i -t j \^dt 1 ---dt n . us) 
J° i=l l<K]<n 



2) In [14] it was discovered by Maple calculations the f^ N satisfy Fuchsian 
differential equation with a factorized "Russian doll" structure 

F 2n f% N = with F 2n = L 2n+1 (N)---L 3 (N)-L 1 (N) (19) 

F 2 „+i/ (2n+1) = with F 2n +i = L 2n +2{N) ■■■L 4 (N)- L 2 {N) (20) 

where Lj{N) are linear differential operators of order j. 

3) It was also discovered in [T3] by Maple calculations that the operators F n 
have in addition a direct sum decomposition 

F 2n = M 2n +i{N) © • • • © M 3 (N) © Mi (AO (21) 
F 2n +i = M 2n+2 {N) © • • • © M 4 (A) © M 2 [N) (22) 

4) Furthermore, the /jy'jvW have a factorization property first found in |14j 
by computer computations and proven for n = 1, 2, 3 in |15j that 

n— 1 2n 
m— m— 
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,(2n+l)/.% n-1 ,(2m+l)/.s 

{jVvjV W _ 7f(2n+l)/Ar\ ^,jV W 

m— 

2n+l 

+ E C , (?" +1 )(7V;i) • • Fff +1 , (24) 

m=0 

where Fn is the hypergeometric function 

F N = 2 F 1 (l/2,N+l/2;N+l;t), (25) 
and /^fjv = 1- The Km\N) depend only on N and we note in particular that 

K^(0) = l (26) 

4 4 \0) = 0, 4 4) (0) = i (27) 

^(0) = -^, tf< 5) (0) = ± (28) 

^ 6) (0) = 0, K[ e \0) = -1, 4 6 \0) = 1 (29) 

The Cm (-^5 i) are polynomials in i of degree for N > 1 

degCj>>(tf;i) = degC r (f +1 >(iV;i) - n-(2iV+l), (30) 

with Cm\N;t) ~ i m as i ~ 0. which have the palindromic property 

C^(iV;t) - H 2A,+1 )+ m - ci 2n) (iV;lA), (31) 
C(2" +1) (iV;0 = t^ 2N +^+ m ■ d 2 n n+1) {N-l/t). (32) 

Explicit formulas for the polynomials Cm' (JV, t) have been obtained in [15] for 
n = 1,2,3 and conjectured for n = 4. For example i-Tp = iV/2 and 



C%\N;t) = (-l) m+1 -! 






' (2N + 


1) 2 1 




4N(N 


+ 1)J 



m 2JV+1— to 



E <$»W> ( 33 ) 

where for < n < N — 1 

^ n (N) = c^ N _ 1 _ n (N) = YZ= cik(N)a n - k (N), (34) 

tfl(N) = c^ 2N _ n (N) = ELo a k (N)a n ^ k (N + 1), (35) 
and for < n < N 

c^n(N) = c^ 2N+1 _ n (N) = ci%(N+l), (36) 
and 2 

#M N ) = (^T ) + 2 ^(l/2)} (37) 
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where 

_ (l/2Ml/2-iV)„ 
an[N) - (1-N) n n\ (38) 

and 

N-l 

^(V2) = E ^ (39) 

It is certainly true (but not yet proven) that the factorizations (|23|) and (|24|) 
hold for all f^ N - The computations in [IS] are based on Fuchsian differential 

equations for the f^\f(t). For n = 4 the order of these equations is 20. These 
equations have a direct sum decomposition into operators which are homomor- 
phic to symmetric powers and products of the operator which annihilates the 
hypergeometric function F^. 

It is furthermore very suggestive that this factorization property has been 
previously seen in the correlation functions of the XXZ model [TB]-|22j. 

The final property of the form factors to be discussed can best be illustrated 
by making a "lambda extension" , first introduced in [23] , of the expansions (jTTJ) 
and (p~2|) by defining 

oo 

C_(M,iV; A) = (1 - tf/^l + £ \ 2n f£%} ( 4 °) 

and 

oo 

C + (M, N ; A) = (1 - t)V4 £ A 2 "/^, (41) 

which reduce to the Ising correlations below and above T c when A = 1. By 
use of a remarkable set of relations presented by Orrick, Nickel, Guttmann and 
Perk [24] in 2001 for small values of M and N, these lambda extensions can be 
written in terms of theta functions |25] 



3 (u; q) = 1 + 2 I"* cos2nu (42) 

n=l 
oo 

e 2 (u; q) = 2q 1 / i q n( - n+ ^ cos[(2n + l)tt] = gVV tt 8 (u + ttt/2; q)(43) 

?i=0 

and their derivatives 

^-e n (u-q) = 6' n (u-q) (44) 

where t x l 2 = k is the modulus of elliptic functions which is related to the nome 
q by 

q = e -«K'(t 1/2 )/K(t 1/2 ) (45) 

and 

i^ 1/2 ) = f 2 f\(l/2,l/2;0 (46) 
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is the complete elliptic integral of the first kind with K'{t 1 / 2 ) = K((l — t) 1 / 2 ) 

The simplest example given in |14j is for the low temperature case with 
M = N = 

C_(0,0;A) = y,"' g ? where A = cosm. (47) 

#3(0; q) 

For the special values A = cos(7rm/n) we find that C_(0, 0; A) and t satisfy an 
algebraic equation. Calling C_(0, 0; A) = r, it is seen in [14] that for A = cos7r/3 

16r 12 - 16r 8 - 8(t - l)r 3 + t(l - i) = 0, (48) 

which is a curve of genus one. For A = cos(7r/4), 

16r 16 + 16(t-l)r 8 + t 2 (t-l) = (49) 

is a curve of genus three which has the simple algebraic expression 

C_(0, 0;cos(^/4)) = 2~ 1 / 4 {1 - t)V"[i + (1 - i) 1 / 2 ] 1 / 4 (50) 

Further results in this direction are [14] 

CL(1, 1;cos(tt/4)) = 2~ 3 / 4 (l - t)V"[i + (1 - i) 1 / 2 ] 3 / 4 (51) 
CL(2, 2;cos(7r/4)) = 2~ 5 / 4 (l - i)V"[i + (i _ t)Va]S/4p _ ^ _ t ) l / 2 }/A 

(52) 

Further results which follow from [24] are given in [26] 

C + (0,0;A,^ (53, 



smu6 2 (0;q)6 2 {0;q) 



C-(W) = - ^"?r*L n ._, (54) 



C + (1,1;A) = - . .y^.. . (55) 

where is to be noted (for N = 0,1) that C+ (N, N; A) is obtained from C_ (JV, TV; A) 
by the interchange 9 2 <-> #3- 

Many further results for various low values of M, TV remain (in the tradition 
of Kaufman and Onsager) to be published by the authors of [14] . 



4 Jimbo, Miwa and Painleve 

The immediate object of the computation of the leading term in the form factor 
expansion by Wu [7] for the row correlation (cro,oCo,Af) and by Cheng and Wu 
[TT] for the general case (t7o,o cr A/,A') was to compute the leading behavior of the 
correlations functions for large separations R = (M 2 + TV 2 ) 1 / 2 . They found that 
for T < T c the correlation decays to the limiting value © as 

(* ,o*m,n) ~ (1 - 1/4 {1 - ^V^- (T) } (56) 
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and vanishes for T > T c as 

(a , a M ,N) ~ (1 - t) 1/4 ^? e -^+ (T) , (57) 

where, in addition to depending on the temperature T, the R independent 
quantities C±(T) and £±(i) depend on the ratio M/N. It is found in [7J and 
[IT] as i ->■ 1 that 

£fc(T) - ^ (58) 



R) 5 



C -( T )-7^W (59) 



and 



where again the amplitudes ^4{.± and A± depend on the ratio M/N. Neither of 
these asymptotic leading terms reduces to the result valid for T = T C (i.e. t = 1) 
where in [7J Wu found that the diagonal correlation has the leading behavior 
for large N 

At 

(<Jo,oVn,n) ~ JjTJI ( 61 ) 

and 

A Ta = 2 1 /i2 e 3C(-D ( 62 ) 

with C'( — 1) the derivative of Riemann's zeta function at —1. 

The history of the result (|6ip is romantic in its own way. In the original 1949 
paper of [5] there is a remark that the diagonal correlation vanishes "slowly" . 
In 1959 Fisher [27} derived the exponent 1/4 and remarked in footnote 8 that 

Onsager, private communication, has derived exact expressions for the cor- 
relations along the main diagonal • • • 

This computation was never published and perhaps there is another typescript 
out there waiting to be discovered. 

Wu [7J also found the large N behavior of the row correlation (oq ) o<7o,.!v)> 
which has the same dependence on N as (|61[) but with an amplitude 



A row = A Tc (cosh 2E h /k B T c ) 1 '\ (63) 

The first purpose of the paper [8] was to connect the three different asymp- 
totic behaviors (|56j) . ([57]) and (f6"Tj) by defining an interpolating function, tradi- 
tionally called a scaling function, 

G±(r)= lim (1 - i) _1/ Vo,o<7Af.;v) (64) 

M.N— >oo,t— >■! 



with 

f anh2E h /k B Tc \ 1/2 ^ f s\xv\v2E v /k B T c \ 1/2 ^ 2 
\smh2E v /k B T c ) \smh2E h /k B T c ' 



! 1/2 

(1 — t) = r fixed. 

(65) 
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For this purpose the exponential representation of the correlation functions was 
derived. When the scaling function was computed it was discovered that G±(r) 
is expressed in terms of a Painleve equation of the third kind 

d 2 r\ 1 { drj\ 2 1 drj 3 



dp-r,\d6j ode' 7 *- 11 ' (66) 

as 

G±( r ) = o i /0 M/2 ex P7 / m Z, ' ( 67 ) 



277^/2)1/2 ex P 4 
with the boundary condition 



r/2 V 



7){e) ~ 1 - -\K (26) as 9 — ^ oo, (68) 



7T 



where Kq(20) is the modified Bcssel function and A = 1. 

This result was first announced in [28 and [29 . Two different proofs were 
given. The first, in [5], is based on Myers' work [3D] on the scattering of elec- 
tromagnetic radiation from a strip and the second |23j is based on a direct 
manipulation of the exponential representation in the scaling limit. 

It is at this point that I first learned of the existence of Sato, Miwa and 
Jimbo when in 1977 I received in the mail (how long ago it was that papers 
were sent by mail) a letter by the three of them with title "Studies on holonomic 
quantum fields II" [31. which generalized several of the results of [5] and made 
clear the relation of the Painleve III equation with the massive Dirac equation. 
This letter was followed by many more where the only change in the title was 
that the Roman numeral was different and by a series of 5 papers with the 
title "Holonomic quantum field theory" [32] . These papers culminated in the 
groundbreaking paper "Studies on holonomic quantum fields XVII" [33) where it 
is derived that the diagonal Ising correlation function for a general temperature 
on the lattice and not in the scaling limit satisfies the sigma form of the Painleve 
VI equation 



N 



The diagonal correlation is related to a for T > T c by 

a(t) =t(t-l)-- \og((To fi a NtN ) - 1/4, (70) 
with the boundary condition at t = of 

Kcw) = t N ^ { -^L + (t^), (71) 



and for T < T c by 



a(t) = t{t - 1) • - \og(a ^N, N ) - t/4 (72) 
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with the boundary condition 

(a , m) = (1 - t)V\l - ^ (^ti) 2 + 0(t^) } (73) 

where (a) at = a(a + 1) ••• (a + iV — 1) for 1 < iV and (a)o = 1 is Pochammer's 
symbol. These boundary conditions are obtained from the leading terms of (|16|) 
and (flT|) as i — >• 0. Furthermore the lambda extensions (|40|) and (|4Tj) satisfy the 
same Painleve VI equation ([69]) where the A appears as a boundary condition. 

The six Painleve equations have a long history [34] , [35] . They are defined as 
those second-order nonlinear equations the location of whose branch points and 
essential singularities (but not poles) are independent of the boundary condi- 
tions and which cannot be reduced to simpler functions. Painleve obtained three 
of these equations |36j and Gambier |37j obtained the remaining three including 
the PVI equation which in the general case has four parameters. However, the 
specific case of Painleve VI needed for the Ising model (j6"9"]) had already been 
obtained by Picard [38! in 1889. Subsequent to the discovery that this PVI 
equation characterizes the diagonal Ising model, this equation has appeared in 
many contexts |39]-|41j ranging from Poncelet polygons to mirror symmetry. 
The sigma form of the Painleve equations was first obtained by Okamoto [42] . 



5 The susceptibility 

The second purpose of the paper [8] was to begin the study of the magnetic 
susceptibility at zero magnetic field X (T), which is computed in terms of the 
correlation functions as 

oo oo 

k B T x (T)= J2 E {(°o,o<jm,n)-M 2 } : (74) 

M=— oo AT=— oo 

where M 2 is the square of the spontaneous magnetization which was given in 
([9]). In order to evaluate the sums in ([74]) the exponential forms (prf and (fT5"j) 
which were the basis of computing the Painleve III equation cannot be used and 
instead the exponentials must be expanded into the form factor representations 
pd|) and (fl"2"|) . Using these forms the sums over M and N are easily evaluated 
as geometric series and the susceptibility is written as the infinite sum of n 
"particle" contributions 

oo 

k B T X +(T) = (l-tf/H-^Y.X^HT) for T > T c (75) 

3=0 

OO 

k B T X -(T) = (l-t) 1/4 5> (!W} Cr) forT<T c . (76) 

In [8] the terms x {n) { T ) for n = 1,2,3,4 were studied. In the scaling limit 
the scaled X (T) for general n were given by Nappi [43] in 1978. For arbitrary 



11 



temperature the results in the isotropic case were obtained by Nickel [5] and 
[TO] and for E v ^ E h in [24] 



X W (T) = 



cot-' a /"^ da>i /' 7r 



with 



where 



2tt 



2tt 



l' n = cot a 



£ — cos w n — — cosa; n ) 2 — (cot a)" 



sinh7 n = cot 2 — cosw„) 2 — (cot a) 4 , 



cot a = Sft/s„ 

e = (i + .s,; 2 ) 1/2 (i + a i/2 

s„ = sinh 2E v /k B T s h = sinh 2E h /k B T 



(77) 
(78) 
(79) 

(80) 
(81) 
(82) 



with 



=[ H h lk 

^l<i<k<j 

sin - wjt) 1 sinh ^(7^ - 7*.) 

/i <fe = cot a , \ = — ■ -, 

sinh 2 \li - Ik ) cot a sm - (tjj + uj k ) 



(83) 



(84) 



and u>j is defined in terms of the remaining Wi from u>i -| ujj =0 mod 2ir. We 

note in particular that for E v = E h 



x (1 \t) = 



t i/i 



with t given by ([T^]) and 

x (2) W 

with i given by ([TU|). 



(1-^/4)2 

(l + t)^ 1 / 2 )-^-*)^ 1 / 2 ) 

3^(l-i 1 /2)(i_ i ) 



(85) 



(86) 



5.1 The amplitude of the susceptibility divergence 

The study of the susceptibility from the form factor expansions was initiated 
in 1973 in 28^ where it was demonstrated that as T — > T c ± the susceptibility 
diverges as 



k B T X (T) ± ~ C ± 
where in the isotropic case 



-7/4 



n/2. 



s = smh2E /k B T. 



(87) 
(88) 
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The constants C_ and C+ arc different and are given as infinite scries 



c - = J2 ° {2n) C + = J2 c(2n+1) ( 89 ) 

n— 1 n— 

where the are n— fold integrals coming from the form factor expansion and 
have been studied both numerically for n = 1, • • • , 5 |28j . [8] . The first term in 
each of (f8"9"]) has been analytically evaluated in [25].[5] 

C« = 1 = -L (90) 

127T 

and the next leading term was evaluated by Tracy [44] as 

C {3) = ^2 (y + 2 - 3^Cl 2 (vr/3)) (91) 

where 

CW») = f;^ (») 

n=l 

is Clausen's function and 

In the tradition of Onsager and Kaufman [3] the details are only in an unpub- 
lished typescript. A curious feature of these results is that the ratio C+/C_ 
is found to be closely approximated by 12tt and the second terms are approx- 
imately three orders of magnitude less than the leading term. The study of 
the constants C_ and C+ has been continued by high precision numerical com- 
putations [23] and the most recent evaluation [5S] in 2011 is to an incredible 
104 places. This is one of the most precisely determined constants in all of 
mathematical physics. 

However, the x i w ) have singularities at other points besides sinh E/ksT = 
±1 and the determination of the analytic properties of the magnetic suscepti- 
bility as a function of temperature has become the most challenging problem is 
the field. 



5.2 Nickel singularities and the natural boundary conjec- 
ture 

The first studies of analytic properties after the initial computations of [S] were 
made in 1999 [5] and 2000 [ID] when Nickel demonstrated for the isotropic case 
E" = E h = E that the integrals ([77]) have singularities in the complex T 
plane on the curve 

|sinh2£/fc s T| = 1, (94) 
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which is the same curve on which the four Pfaffians of Kaufman's original evalu- 
ation [4] of the Ising partition function vanish. This was extended to the general 
case E v 7^ E h in [21] where the singularities of \ {T) are at 

cosh 2E v /k B Tcosh 2E h /k B T 

~smh2E h /k B Tcos(2irj/n) - sinh2E v /k B T cos(2irk/n) = 0. (95) 

with 

< j, k < [n/2], j = k = excluded (96) 

where [x] is the integer part of x and for n even j + k = n/2 is also excluded. 
In terms of the variable used in [46]- [55] for the isotropic lattice with s given by 

dSHD 

w' 1 = 2{s + S - 1 ) (97) 
these singularities for n = 3, 4, 5, 6 are given in table 15.21 

Table 1: The Nickel singularities of x^ f° r n — 3,4,5,6 



n 


w 


3 


-1/2, 1 


4 


±1/2 


5 


-i -1±V5 3±v / 5 
> 4 > 2 


6 


±1, ±1/3 



where we note that sinh 2E/k B T is real for — 1/4 < w < 1/4 and is complex with 
| sinh2_E/fcBT| = 1 for 1/4 < |io|. If we call e the deviation from the singular 

temperatures T^z , determined by (J95J) , then for T > T c the singularity in 

x (2j + l)( T ) is 

e 2 j(j +i)-i lne (98) 

and for T < T c the singularity in ^^(T) is 

e 2 ^ 2 - 3 / 2 . (99) 

It is striking that the number of singularities increases with n and becomes 
dense in the limit n —> oo. This feature led Nickel to the conclusion that unless 
cancellations occur there will be a natural boundary in the susceptibility x(T) in 
the complex T plane at the location (|95l) . The existence of a natural boundary 
in the complex temperature plane is not contemplated in the scaling theory of 
critical phenomena. 

5.3 Puchsian equations 

The next step in the study of the susceptibility was begun in 2005 [46] and has 
continued in the series of papers [47] -[55]. In these papers exact Fuchsian differ- 
ential equations for the X (T) in the isotropic case E v = E h are determined 
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by use of Maple by first expanding the integrals in an appropriate variable such 
as w or w 2 and then using Maple programs which obtain ODE's from these 
series. The resulting differential equations have very special properties such as 
being globally nilpotent [32] which allow for extensive analysis to be carried 
out. These studies have uncovered several new and important features of the 
susceptibility; namely that the X^ n \ w ) have a direct sum decomposition and 
that they have further singularities beyond those of (|95|) . 



5.3.1 Direct sum decompositions 

In [2PJ and [52] it is shown for 1 < n < 6 that x^ ( w ) have the same direct sum 
decomposition seen already in the diagonal form factors 

X (2n) (w) = ]T K^ X {2m) (w) + rt 2n \w) (100) 

m— 1 

n-l 

x (2n+1) (w) = Y, ^ 2n+1) x (2m+1) M + ^ (2n+1) M (ioi) 

m— 1 

where the ^ n \w) satisfy Fuchsian equations of order m 

L (n) . jj(n) = o (102) 

with 



n 3 4 5 6 
m 6 8 29 46 



(103) 



(n) 

The Kj are constants which for n — 3,4,5,6 coincide with the values of 
K£\0) given in f26>f29]). 

The operators in (|102p factorize further. For Lg 3 -* and L^f 1 we have we have 

4 3) = 4 3) • 4 3) ■ 4 3) (104) 

and 

4 4) = 4 4 M 4 ^4>4>4S) d05) 

where the numeral in the subscript indicates the order of the operators which 
are given in [25] and [25]. The operator L^g has been found in [5T],[53] and [55] 
to have the factorization 

L$ = Lf ■ ■ I? ■ 2# (106) 
where 4i^ has the further direct sum decomposition (A.l) of |55] 

4? = (Z 2 ■ Ni) ®V 2 ® (F 3 • F 2 ■ L{) (107) 
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Similarly in (56) and (57) of [54] the operator L^ 6 ' is shown to have the decom- 
position 

I^=Lf.L$.L$ (108) 

where has a direct sum decomposition into the sum of four operators but the 
possible reducibility of has not yet been determined due to computational 
complexity. 



5.3.2 Singularities 

The location of the singularities of the operators are obtained by examining 
the roots of the polynomial multiplying the highest derivative d m /dw m and this 
analysis shows that there are further singularities beyond the singularities at 
w = ±1/4, oo and the Nickel singularities (|95]) . 

In [47] that the differential equation for x^ 3 ' ( w ) admits additional singular- 
ities at 

w = t-^— (109) 



which correspond to 



-l±iy/7 



s\ = V2 (111) 

where we note the the singularity at (|110j) is inside the unit circle |s|=l and thus 
cannot appear in the principle sheet of the integral for x^ which is analytic for 

N<i. 

There are no additional singularities in x^ ( w ) and the singularities of 
X^ 5 \w) are shown in (34) of [51] to be at the roots of following polynomial 



w 



(1 - 4w) 22 (l + 4w) 16 (l - w) 4 (l + 2w) 4 (l + 3w + Aw 2 ) 



(1 + w)(l -3w + w 2 )(l + 2w - Aw 2 ) 

(1 - w - 3w 2 + Aw 3 )(l + 8w + 20w 2 + 15w 3 + 4w 4 )(l - 7w + 5w 2 - Aw 3 ) 
(1 + Aw + 8w 2 )(l - 2w) (112) 

The singularities located by the roots of the first line in (|112l) are identical with 
the location of singularities of x^ and the roots of the second line are the 
Nickel singularities of of X ■ Most of the remaining singularities correspond to 
complex values of s not on |s| = 1. 



6 Diagonal susceptibility 

The integrals (|77|) for the n particle contribution to the susceptibility X (T) are 
quite complex and the Maple-based studies cannot be extended much beyond 
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their present limits. Therefore it would be of great utility if a simpler set of 
integrals could be found which would still incorporate all significant analytic 
features of the the x^ ■ Several such simplified modifications of the integrals 
have been studied [5U] but by far the most natural case is to restrict the two 
dimensional sum over the lattice positions M, N in (|74"|) to the lattice diagonal 
M = N and thus to consider the susceptibility that will result if a magnetic 
field is applied only to the diagonal 

oo 

k B T X d{t)= {(vo,o<tnm)~M 2 } 7 (113) 

N=—oo 

where the dependence on T is now for all E v and E h in terms of the single 
variable t defined by (TO} for T < T c and by (TO} for T > T r . 

This diagonal susceptibility has been studied in [55] and [57] and has been 
found to have the remarkable simplification over the bulk susceptibility that all 
singularities of the differential equations are at s — 0, oo and \s\ = 1. There are 
no other complex singularities for \s\ ^ 1 such as appear in X ^(t)- Further- 
more x^ 3 -* (t) and Xd (*) have been found to be explicitly expressed in terms of 
generalized hypergeometric functions p +iF p . 

6.1 Integral representations 

From the integral expressions for f^ N (t) of (TO} and (T7} given in [TO] and [TO] , 
we find in [53] the expansion for T < T c 



and for T > T c 



oo 

k B T X dAt) = (1 - t) m E xf n+1) (t), (115) 

n=0 



where 



xfif v /2 n n (i-*^-i*«)- a 

J-J- V«2i(l-a;2j-i)(l -tx 2j -i)J }i ,Ai 

j=l v J l<j<n l<fc<n 

x JJ (#2.7-1 - 2J2fc-i) 2 (x 23 - a; 2fe ) 2 (116) 

l<j<fe<n 



and for T > T c 



fe=i 
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1 +t n + 1 / 2 x 1 ---x 2n+1 fw Ni/a 



-1/2 



1 - - /'<+' 2 .r> ■■■ • - 

J'=l 

n+1 

x II (j 1 _a; 2j-i)(l -ia^-i)-^-^ 
i=i 

x n n ( x -^2j-n2*)" 2 

l<i<n+l l<fe<n 

X fj (%2j-l ~ X 2 k-l) 2 Yl ( X 2j-X2k) 2 - (H7) 

l<j<k<n+l l<j<k<n 

The expressions (|116|l and fll 1 7[) are, indeed, much simpler than the correspond- 
ing expressions for x^ given in (|77l) . In particular 

and 

which are simpler than (1551) and (|86l) respectively. Most noticeable is that 
x' 2 )(w) in (|86p has a logarithmic singularity at i = 1 (ic = 1/4) while Xa (*) m 
(fTTg|) does not. 

6.2 Root of unity singularities 

In addition to the singularity at t = 1 it is straightforward to see from the 
integral expressions (111 6[) and (|117l) that Xd*^ (f) nas singularities at 

to = 1 (120) 

of the form 



,2n -1 



€ 



hie (121) 



and Xd* n+1 \t) nas singularities 

t n +1/2 = 1 (122) 

of the form 

£ («+l) 2 -l/2 (123) 

where e is the deviation from to- These are the analogues for the diagonal 
susceptibility of the Nickel singularities of the bulk susceptibility x of (|95|) . 
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6.3 Direct sum decomposition 

The x [ d ] (*) have the 

same direct sum decomposition seen already in the diagonal 
form factors and x ( w ) 

n-l 

X ^\t) = £ Kf.- ] xf\t) + ti d 2n \t) (124) 
i=i 

xT +1 \t) = £ ^ +1) xf +l) {i) + nr +1) (t) (125) 

where ^Q™ are constants. However, unlike X^( w )i the operators L d n ^ which 
annihilate (t) have a further direct sum decomposition 

L ( &=L$ + L$ and L$ = L% + L$ (126) 

(3), 



6.4 Results for x«i (*) 

For xf^ we explicitly hnd by combining [5 2) and |56j and setting a; = f 1//2 
that 

where 

X 

I .. . .. .. , o, I 

X d -M x ) = 7T~ 

and 



*f (*) = sxS(s) + o^S(^) " «X$l(x) (127) 



xg(-) = T ^=xi 1) W (128) 
xg(z) - TfT^aJiClA -1/2; 1; x 2 ) - ^^(1/2, 1/2; 1; .t 2 ) (129) 



20 

+ -^F(l/6, 1/3; 1; Q)(F(7/6, 4/3; 2; Q) 



(130) 



with 



^ " 4 (a: 2 +x + l)3 



where we note that 



(1-^(1 + 2^(2+^ 
1 Q ~ 4(l + x + 2 ; 2 )3 

From (fHT)) and (fl"2"g|) we see that in (fl"2"5"|) we have i^ 1 ' = 1/3. 

We see from (|117[) that xi 3 ^( x ) vanishes when x — > as a; 4 . However, ^^(a:) 

(3) (3) 

and X^.^x) are constant as x — > and Xd-2( x ) vanishes linearly in x. The three 
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constants in (jl2T[) are determined by matching with the x A behavior of x^ and 
this requires that the three constants will solve a set of five (overdetermined) 
linear equations. 



As x — > 1 we find that diverges as 



(3), x 1 (\ 5tt 4tt \ 

Xd {> i- x \3 i8r 2 (5/6)r 2 (2/3) r 2 (f/6)r 2 (i/3)y 



0.016329- •• 



1-x 

(133) 

Furthermore X&a{ x ) nas an additional singularity at x — > e ±27 ™/ 3 which, to 



!;3 

leading order is 



6.5 Results for Xd^(t) 



These results have been extended in [56] and [57] to Xd (*) wnere is is shown 
that 

X<£\t) = l X gW + 3^5xg(*) " ^xgW (135) 

where 

xg(*)=X d 2) (*) (136) 

xg(*) = ^p^l/2, -1/2; 1; tf - 2 F X (1/2, 1/2; 1, tf 
2t 

F(l/2, 1/2; 1; t) 2 *i(l/2, -1/2; 1; t) (137) 



and 
with 



1-t 



= ^3 • 4^(11/2, 1/2, 1/2, 1/2]; [1, 1, l}t 2 ) (138) 



2 16i 2 - t - 11 

3 <- 1 
1 31t 2 -At -11 



A 3 = 2(1 + t)t 3 Df + — t 2 d 2 t 



tD t + t (139) 



3 t - 1 

The singular behavior as t — > 1 of Xct\(t) arL d X&a ^ s easily obtained and 

the singularity of x&i at i = 1 is obtained by use of the analytic continuation 
formula of Biihring [6T] . The final result [57] is that as t —> 1 

7 mA-T^ln 2 A d40) 



16tt 2 1 - 1 16tt 2 1 - 1 
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where 

3Ii - 4J 2 = -2.2128121 •• • (141) 

has been given to 100 digits. 

At the root of unit singularity t = — 1 the leading singular behavior is 

268*^ + *) 7 M1 + *) ( 142 ) 

6-6 

The ODE satisfied by x^ ( x ) nas been studied in [57] modulo a large prime. It 

(5) 

is found that the minimal order ODE is of order 19 and that the operator L d .[ g 
has the decomposition 

L d % = L% 8 L% (143) 

where L d ^) 2 is the second order operator which annihilates x^d-2 ( x ) ano - ^d\r nas 
singularities at x = 0, oo, 1, — 1, £3 = e ± 2 W3 ; X5 = e ±27n/5^ e ±47ri/5 w ]j ere 
the non-integer exponents at X3 are 5/2,7/2,7/2 and at £5 are 23/2. It has 
been further found that 

T (5) _ r(5) . r(5) z,^ 

with 



r(5) _ r (3) 



t 

are all given in [5 7) . 



(3) • • (3) 

where L d . m annihilates Xd-m ano - ^ ne remaining operators in this decomposition 



6.7 Singularities and cancellations 

By examining the integral representations for x^ n \ w ) an d XT\t) <|H6[) . 
(| 1 1 7[) it is clear that these integrals have no singularities for | s'mh2E /ksT\ < 
1. The singularities at | sinh2_E/fcsT| < 1 of the differential equations for 
^(«) -will only appear in analytic continuations of the integral in the complex 
plane of the variable w. The corresponding differential equations for X d ft) are 
significantly simpler because they have singularities only at t — 0, 00 and \t\ = 1. 

It remains to discuss the singularities in the differential equations which do 
lie on I sinh 2£ , /fcsT| — 1 and to give an explanation for the observation that the 
singularities of the ODEs for ^("- 2m ) ( w ) an d ^j™ 2m '(£) are also singularities 
°f X^ n \ w ) an d X ft) respectively even though the integrands are singular only 
at the points given by ([95]) for X (n) M an d by (fT20]) and ([T22]) for x<j n) (t). 

The resolution of this is easily seen for j;v(t). By an examination of the 
integrals (|116[) and (|117p we see that there are paths of analytic continuation 
possible in the complex t plane where the contour of integration must be de- 
formed past the pole at 

1 - tx 2 jx 2 k+i = (146) 
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and the residue at that pole will reduce the denominators in and Xd"^ W 

from 

1 - t n xi ■ ■ ■ x 2n (147) 

and 

l-t n+l / 2 Xf-x 2n+1 (148) 

to the denominators in x^ n 2 ^ (t) and Xd ^ (t) respectively with n — >■ n— 1 and 
two less integration variables. Therefore, the singularities of \d W wu ^ n °t 
appear on the principle sheet of the integral which is analytic at t = but only 
on analytic continuations to non-physical branches. The similar phenomenon 
occurs for x^ n \ w )- 

It remains to reconcile this non appearance of the singularities of x^ 2 W 
in the physical sheet of Xd^W with the direct sum decompositions (I124p and 
(|125[) . This will be accomplished by showing that the term fifj (t) has sin- 
gularities which exactly cancel the singularities on This requires the 
solution of a global connection problem which has not yet been explicitly done 
even though from the examination of the original integral the resulting exact 
cancellation must hold. 



7 Conclusion 

Now that we have summarized the known features of the Ising correlations, form 
factors and susceptibility we can proceed to discuss what is not known. This is 
the fascinating, mysterious and thus romantic part of the subject. 

7.1 Conformal and quantum field theory 

One of the most important features of the Ising model is that the scaling limit 
satisfies all the axioms for a massive Euclidean quantum field theory and that 
at T = T c the long range correlations are those of a conformal field theory with 
central charge c = 1/2. This is in fact the earliest conformal field theory known 
and from this beginning a vast new field of mathematics and physics has been 
developed in the last 30 years. However, the Ising model is much more than 
a conformal field theory because we have a vast number of results for T ^ T c 
which are the simplest examples of properties of massive Euclidean quantum 
field theories. Part of the romance is the exploration of how these Ising results 
can be used to extend massless conformal field theories into the massive region. 

7.2 Form factors, exponential forms and amplitudes 

The derivation [12] of the exponential and form factor expansion for the diago- 
nal Ising correlation is much more general than this special case. Indeed in [12] 
it is proven that every Toeplitz determinant Q with a generating function <p(£) 
such that In </>(£) is continuous and periodic on |£| = 1 has both an exponential 
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and a form factor expansion. Furthermore these Toeplitz determinants are also 
expressible as Fredholm determinants [58] (at times in several different ways 
[55]). Consequently the Ising computations have subsequently been extended 
to several very important problems including the seminal work on the one di- 
mensional impenetrable Bose gas and on random matrices by Jimbo, Miwa, 
Mori and Sato [50]. 

To illustrate the differences between the form factor and the exponential rep- 
resentation of the correlation functions, we consider the computation by Tracy 
[52"] of the constant At c of ([52)1 . In the scaling limit the scaled correlations 
in the general case where E v ^ E h depend only on the single variable r (|55]) . 
Therefore we can restrict attention to the scaled form of the diagonal correla- 
tion ((TOjOCjv.jv) an d consider the lambda extension of the scaling form of the 
exponential form (|I4I) which we write as 

oo 

G_(r;A)-cxp^A 2 'V 2n) M I (149) 

n=l 

where 

g( 2n \r) = lim F$>{t), (150) 

scaling 

which depends on the single variable r instead of the two independent variables 
N and t. Tracy finds that, as r — >• 0, 

g^(r) = -a n \nr + l3 n + o(l). (151) 

Therefore, defining the lambda dependent sums 

oo oo 

a(A) = ^ 2 "«« m = E x2n ^ ( 152 ) 

n— 1 n— 1 

we find 

e /9(A) 

G_(r; A) ~ exp{— a(A) lnr + /3(A)} = (153) 

In the Ising case where A = I the functions specialize to 

a(l) = l/4 P(l) = 1nA (154) 

where A is the constant in (|f)T]) . 

If, however, instead of the scaled exponential form we define the scaled limit 
of the form factors /tvatW as 

/( 2 ")(r) = lim /<$(*), (155) 

scaling 

then as r —> 

n 

/( 2 ™)( r ) = 4 2l ° lnfe r + o(l). (156) 

fc=0 
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Thus, in order for f|153[) to agree with the r — > behavior of the form factor 
expansion, we need 



MX) 



oo - / oo \ k oo 1 / oo \ ' 

fc=0 \ n=l / A;=0 \n=l / 
oo n 

1 + E A2n E a i 2n) ln " r 



(157) 



to hold term by term for each power A 2 ". The requires an infinite number of 
identities between the a^ n ^ . 

As an additional remark we note that if we rewrite the integral (|16[) for ^ 
as a contour integral, rescale the variables Xk by Xk — t~ x / 2 yk and then send 
U2k —> l/U2k we see that as t — >• 1 the integral has logarithmic divergences as 
in (|156[) . The amplitudes a n are closely related to the special case with p = 1 
of the integral found by Dotsenko and Fateev [53] in their study of four point 
correlations in conformal field theories with central charge c < 1 



I n ,m(a, P; p) 



X 



n\m\ y Q jla y Q 

n^-*i)i v in^-^)i 2p ii 



where P indicates the principal value and 

a' = -p'a, P' = -p'l3, p' 



(158) 



(159) 



7.3 Exponentiation 

Form factor expansions exist for many massive models of quantum field theory 
including sine-Gordon and the non-linear sigma model [64] and similar form 
factor expansions exist [55], [55] for the XXZ model on a chain of finite length 

L 

Hxxz = E + + Aff K+i + H °j} ( 16 °) 

i=i 

Moreover the Feynman expansion of amplitudes in quantum field theory is also 
what we have called here a form factor expansion. In all of the models there 
are limiting cases where series of multiple dimensional integrals expand to se- 
ries in powers of logarithms which need to be summed. However, unlike the 
Ising correlation functions these form factor expansions do not come from ei- 
ther Toeplitz or Fredholm determinants and thus the exponentiation methods 
of Ising correlations are not applicable. 

Over the years an immense effort has been made to sum the form factor 
series of logarithms. In quantum field theory this starts with the classic 1939 
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paper of Bloch and Nordsiek [57] on resummation of infrared divergences in 
quantum electrodynamics. A second example is the Regge theory of the 60's 
and 70's an 2nth order Feynman diagram expansion of a four point scattering 
amplitude is shown to diverge as the energy s — ¥ oo with a fixed momentum 
transfer t as 

9 2n a(t) nl ^ (161) 

This is a "leading" log approximation and is analagous to the Ising case if only 
the first term in the series (|152[) for a n is retained. More recently there has been 
a great deal of work on quantum chromodynamics |69| where many non leading 
terms summed by the use of an ingenious decomposition of the multidimensional 
integrals. 

In the theory of integrable systems a great deal of effort has been devoted 
to compute the long range asymptotic behavior of the correlations of the XXZ 
model in the massless region — 1 < A < 1 from multiple integral representations. 
One method is presented in [70] which shows how to modify the Fredholm 
determinant form which holds for A = by suitably picking out the important 
pieces of the multiple integrals. This has led to the computation of both the 
exponents and the amplitude of the long range behavior of the correlations when 
the H ^ 0. The study of the correlations from the form factors is begun in [65] - 
[55] with more results announced to be forth coming. A full exploration of the 
relation of these subjects is beyond the scope of this article. 

7.4 Short distance versus scaling terms 

In the n-particle expansions of the full ([76)) , (|75|) and the diagonal (jl 14[) , (1115j) 

susceptibility the x^ n \t) an d the Xd^if) wn ^ (l° r n — 3) have terms which 
contain powers of Int. From this it might be inferred that the susceptibility will 
contain terms of the form (1 — t) 1 / 4 +P ln 9 (l — t). However, from the extensive 
calculations on long low and high temperature series expansions made in |24j and 
[35] such terms do not appear to exist. Instead the susceptibility is conjectured 
to have the form for t — > 1 of 

[Vol 

k B T X (t)± = (i - r 7/4 e c± (i - ty + E E b ± q) (! - ^ lnP ( x - f ) ( 162 ) 

j=0 9=0 p=0 

The first term is called the "scaling function". The second term is called the 
"background" or "short distance" term and is numerically obtained by summing 
correlation functions instead of form factors. In [45] it is stated that the "scaling 
function" is determined by conformal field theory while for the "short distance" 
term here is "no explicit prediction" . In [23] the belief is stated that the sepa- 
ration into "scaling" and "short distance" parts is "tantamount to the scaling 
argument that in the critical region there is a single length scale proportional to 
(1 — t)~ u with v = 1". It would be highly desirable if this distinction between 
"scaling" and "short distance" terms could be made precise and if both terms 
could be obtained by use of the form factor expansion alone. 



25 



7.5 Natural boundaries and A extensions 

Perhaps the most perplexing question concerning the relation of the Ising model 
on a lattice with the scaling field theory limit is the existence of the natural 
boundary in the susceptibility implied by the singularities ([55]). found by 
Nickel [5], [TU]. The magnetic susceptibility is the second derivative of the free 
energy with respect to an external magnetic field H interacting with the spins 
as —H^2j k a J- k - ^ n me scaling limit the Ising model in a magnetic field is 
also a field theory and the analyticity properties of this field theory have been 
extensively studied by Fonseca and Zamolodchikov [71 with the conclusion that 
there is no natural boundary. How can this be reconciled with the computations 
of [9] and [10]? 

The existence of the natural boundary suggested by Nickel in [9] and [10] 
rests on the accumulation of the singularities (|9"5|) and ([9"9"]) and the assumption 
that there is no cancellation. However, for this argument to hold we need to 
be able to show that the limit of t approaching the location of the supposed 
natural boundary (I95[) will commute with the infinite sum over the n particle 
contributions x (n) ( T ) in CZ3 and (J7BJ). Since the natural boundary does not 
exist if only a finite number of the \ W(T) are included this interchange need 
to be investigated. It is also possible that the existence of a natural boundary 
could depend on the value of A in the lambda extensions of (|75l) and (I7B1 

oo 

k B T X +{T;X) = (l-t) 1 ^- 1 /^^))^ for T > T c (163) 

j=o 

OO 

k B T X ~(T;\) = (l-t) 1/4 5> v X (!W) Cr) for T <T C . (164) 
These possibilities remain to be investigated. 



7.6 Row correlations 

All of the results obtained for the diagonal correlation, which depend on the 
single variable t, can be extended to the row correlation, which depends on the 
two variables a± and a 2 in a symmetric fashion ([5]). In particular it has been 
pointed out to me by Jean-Marie Maillard and Nicholas Witte in private con- 
versations that the Painleve VI results of Jimbo and Miwa [33J can be extended 
to a two variable Gamier system |72j . However, this system must possess some 
most interesting properties because one of the most important properties of the 
Ising model is the fact that, when these two variables are rewritten as 

sinti 9 7^ v Ik 

k = smh2E v /kBTswh2E h /k B T and r = — r 1 ^—. (165) 

sinh2i/' l /KBT 

the dependence on k (the modulus of the elliptic functions) and the anisotropy 
ratio r which is related to the spectral variable of the star triangle equation [73] 
is dramatically different. These results for Gamier systems have also yet to be 
obtained. 



26 



8 Romance versus Understanding 



In a lecture given in Melbourne in January 2006 [74], I gave the following defi- 
nition of "understanding" 

No one can be said to understand a paper unless he is able to generalize the 
paper. 

This definition is open to criticism on at least two grounds. Firstly the use of 
the word "he" has a sexist implication which is neither appropriate nor intended. 
Secondly, there are surely subjects which are fully understood where further gen- 
eralization is pointless. An illustration of this are the laws of thermodynamics 
which have been fully understood by physicists for many decades (even if they 
are not accepted by the overwhelming majority of voters and politicians) . 

However, precisely because thermodynamics is fully understood, it has lost 
the mystery it had at the time of Gibbs, Boltzmann and Ehrenfest. This illus- 
trates the great truth that understanding is the enemy of romance because once 
the mysteries are understood the romance dies. 

Fortunately for romance, there are many mysteries of the Ising model which 
are far from being understood. The romantic in me says that, even when these 
mysteries have been understood, the understanding of the mysteries will gener- 
ate new mysteries and the romance of the Ising model will be everlasting. 
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